Multipartite W states for chains of atoms conveyed through an optical cavity 
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We propose and work out a scheme to generate the entangled W states for a chain of N four-level 
atoms which are transported through an optical cavity by means of an optical lattice. This scheme 
is based on the combined laser-cavity mediated interaction between distant and equally separated 
atoms and works in a completely deterministic way for qubits encoded by two hyperfine levels of the 
atoms. Only two parameters, namely the distance between the atoms and the velocity of the chain, 
determine the effective interaction among the atoms and, therefore, the degree of entanglement that 
is obtained for the overall chain of N qubits. In particular, we work out the parameter regions 
for which the Wat states are generated most reliably for chains of A/" = 2, 3, 4 and 5 atoms. In 
addition, we analyze the sensitivity in the formation of entanglement for such chains of qubits due 
to uncertainties produced by the oscillations of atoms in optical lattices. 

PACS numbers: 42.50.Pq, 42.50.Dv, 03.67.Bg 



INTRODUCTION 



During the last decades, quantum entanglement has 
been found essential not only in studying the non- 
classical behavior of composite systems but also as im- 
portant resource for the engineering and processing of 
quantum information. Nowadays, numerous applications 
are known that greatly benefit from having entangled 
quantum states available, such as super-dense coding [l| , 
quantum cryptography [25, or Grover's quantum search 
algorithm [3] to name just a few. Despite of the re- 
cent progress in dealing with composite quantum sys- 
tems, however, the controlled manipulation of these sys- 
tem and their interaction with the environment remains 
still a great challenge. Beside of various other imple- 
mentations of composite quantum systems for quantum 
control and applications in quantum information, a ex- 
cellent control in generating entangled states has been 
achieved recently with neutral atoms that are coupled to 
a high- finesse optical cavity 0, [11 . 

In practice, there are two typical ways to encode a sin- 
gle qubit into the level structure of an atom: Apart from 
(i) selecting two levels separated by an optical transition 
frequency {optical qubit), one may also utilize two hyper- 
fine levels of-typically-the atomic ground state, some- 
times referred to as hyperfine qubit. In contrast to the 
optical qubits, the use of hyperfine qubits has the ad- 
vantage of long coherence times (~ 1 s) and, moreover, 
such qubits are known to be more robust with regard to 
external perturbations or stray fields. Finally, a num- 
ber of microwave techniques have been developed during 
the last decades that allow to initialize, manipulate, and 
detect the states of such qubits 
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FIG. 1: (Color online) Schematic view of atoms in an op- 
tical lattice (conveyor belt). Two focussed and counter- 
propagating laser beams with frequencies uu-\-A/2 and cj — A/2 
give rise to an interference pattern in the field strength with 
a series of equidistant potential wells in which the atoms can 
be trapped. The distance d between two neighbored wells is 
given by (half of) the lattice wavelength A, while the velocity 
of the belt v is determined by the detuning A of the two laser 
beams. 



However, the hyperfine qubit (s) cannot couple directly 
to a cavity with a resonant mode-frequency in the opti- 
cal domain. Therefore, a four-level configuration need to 
be considered, in which the two hyperfine levels are sup- 
plemented by two electronically excited levels which are 
separated by the optical transition frequency and being 
compatible with the resonant frequency of the cavity. In 
order to realize an effective manipulation of the hyper- 
fine qubit by means of the optical cavity, its state must be 
mapped coherently upon the electronically excited states 
and back to the hyperfine levels, once all the desired 
atom-cavity interactions have been performed. Indeed, 
such a indirect coupling between the hyperfine qubit and 
optical cavity opens a route for the generation of entan- 
glement and more complex quantum states between two 
or more atoms by means of the cavity-mediated interac- 
tions. 



Despite of the recent progress to couple one single atom 
to the cavity mode, further control of the atomic motion 
is necessary in order to manipulate the interaction be- 
tween cavity and a chain of atoms. At the same time, an 
excellent control of the motion of atomic chains is merely 
possible by using optical lattices (conveyor belts) [9|, 
which have recently been utilized in various setups of cav- 
ity QED [10-12]. In Fig. [T] we displayed a schematic view 
of such optical lattice in which two counter-propagating 
laser beams with parallel linear polarization produce an 
interference pattern in the field strength that gives rise 
to a series of equidistant potential wells, where neutral 
atoms can be trapped. These wells allow to control the 
position of atoms with a sub-micrometer precision over 
millimeter distances due to their tight confinement along 
the lattice axis. 

The combination of such a lattice with the (optical) 
cavity QED setup, however, make it necessary to revise 
the evolution of the atom-cavity interaction for a chain of 
atoms that is conveyed by such a lattice through the cav- 
ity. In particular, one need to analyze how the spacing 
between the atoms and velocity of the atomic chain (lat- 
tice) will affect the formation of entangled states between 
the hyperfine qubits. By this revised evolution, moreover, 
the small sample approximation, i.e., when the spacing 
between the atoms is considered negligible in compari- 
son to the cavity waist, should be abandoned and the 
position-dependent effects should be taken into account. 

For two hyperfine qubits (in a four-level configura- 
tion) , we have recently proposed a scheme to generate the 
maximally entangled state ^ (e*'^|0i,l2) + |li,02)) by 
means of the combined laser and cavity mediated inter- 
action [13]. In this reference, we considered the position- 
dependent coupling between the atoms and cavity which 
allowed us to describe the formation of entanglement be- 
tween the two atoms being separated by a macroscopic 
distance. An effective interaction between the hyper- 
fine qubits was achieved if both, the laser and cavity 
fields, are detuned with regard to the atomic transition 
frequencies (see below). In particular, we demonstrated 
explicitly how the degree of entanglement depends on the 
atomic velocities and the spacing between two atoms. 

In the present work, we extend this analysis and pro- 
pose a scheme to generate the entangled W state [15] 
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(e^^|li,02,...,0iv) + ... + |0i,02,...,liv)) (1) 



between the hyperfine qubits of N (four-level) atoms 
that are equally distanced from each other and conveyed 
through the cavity by means of an optical lattice. Simi- 
larly to Ref. [13], this scheme works in a completely de- 
terministic way and is based on the position-dependent 
interaction between distanced atoms which is mediated 
by the cavity and laser fields. The two parameters that 
control this atom-cavity-laser interaction are (i) the ve- 
locity of the atomic chain along the axis of the lattice 
and (ii) the distance between the atoms. For the chains 
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FIG. 2: (Color online) (a) Schematic setup of the experiment. 
A chain of N neutral atoms passes through a pair of Raman 
lasers Li, an optical cavity C with a laser beam L, and a 
second pair of Raman lasers L2. The atoms are supposed to 
move in a chain along the z-axis with a constant velocity v 
such that the chain crosses the cavity at the antinode. Apart 
from the cavity waist that is just one half of the minimum 
width of the cavity radiation field, the cavity is also charac- 
terized by its position-dependent coupling strength g{z). (b) 
The atomic four-level A-type configuration in the Schrodinger 
picture and (c) in the interaction picture. 



consisting of N = 2,3,4 and 5 atoms, we determine the 
velocities and distances for which the initially uncorre- 
lated atoms produce the Wat states most reliably. Apart 
from generation of the W states, we discuss also how the 
proposed scheme can be implemented most efficiently and 
analyze how robust are the entangled states with respect 
to small oscillations in the atomic motion as caused by 
the thermal motion of atoms in the optical lattice. 

The paper is organized as follows. In the next sec- 
tion, we first outline our scheme to entangle the hyper- 
fine qubits of N initially uncorrelated atoms. In Section 
III, we then explain and discuss the effective Hamiltonian 
which describes the atomic evolution; we analyze in par- 
ticular the parameter (regions) in Subsections III.A-D 
for which the states W2, W3, W4, and W5 are generated 
most reliably. In Section IV, we later discuss a few issues 
related to the implementation of our scheme and how it 
is influenced by small oscillations in the motion of the in- 
dividual atoms, while a short summary and outlook are 
given in Section V. 
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II. GENERATION OF THE W STATES 

We shall first explain the basic idea of the proposed 
scheme for generating multipartite W entangled states 
within chains of neutral atoms without going much into 
details. We assume that the atoms are initially in a 
product state and that they are inserted into an optical 
lattice being equally separated by a distance d as dis- 
played in Fig.[2ja). Moreover, the atoms are supposed to 
move with a constant velocity v along the (lattice) 2;-axis 
such that their position vectors fi{t) = {0,0, -\- vt} 
cross the cavity at the anti-node and where denote 
the initial position of the i—th atom. As briefly outlined 
above [cf. Fig. [l], this velocity v and inter-atomic dis- 
tance d can be controlled experimentally by adjusting the 
shift in the frequencies of the two counter-propagating 
laser beams and by selecting a proper wavelength of the 
optical lattice, respectively [9]. 

Each of the N identical atoms represents a (hyperfine) 
qubit in a A- type configuration as displayed in Fig.[2fb), 
in which the two hyperfine states |0) and |1) are supple- 
mented by the electronically excited states |e) and \a) in 
such a way, that the transitions \a) |1) and \a) |0) 
are (electric-dipole) forbidden due to the angular momen- 
tum and parity selectrion rules. Initially, the atoms are 
prepared in the product state 

|li,02,...,07v) = 111) X IO2) X ... X lOiv), (2) 

where the numbering corresponds to the (increasing) co- 
ordinates ^5^, ^2 , . . . , of atoms along the z-axis. By 
this assumption, therefore, the qubits are loaded to the 
cavity in the reverse order, i.e. qubit No. 1 corresponds 
to the last atom inside the chain. 

Just before each atom enters the cavity, its electronic 
population is transferred from the state |0) to the state 

I a) with a pair of (slightly) off-resonant laser beams that 

1 2 

are coupled to the atomic transitions |0) \e) \a). 
This population transfer is known as the two-photon Ra- 
man process and can be implemented, for example, by 
means of two phase- locked laser diodes [6[. Below, we 
shall briefly refer to this transfer in the population as 
the Raman pulse and shall distinguish between the Ra- 
man zones Li and L2 in front and behind the cavity [see 
Fig. Efa)]. In this notation, the Raman pulse L2 is uti- 

2 1 

lized to perform the population transfer \a) \e) |0) 
back to the hyperfine level. Inside the cavity, therefore, 
the product state (|2]) becomes |li, a2, . . . , ^at) since the 
last atom is unaffected by the Raman pulse Li. 

To explain the mechanism of the cavity-mediated in- 
teraction between the atoms in more detail, let us first 
consider a chain of two atoms prepared in the prod- 



uct state |ei,a2), and where each atom is coupled to 
a detuned optical cavity by the transition \a) ^ \e) 
[cf. Fig. EJb)]. In this case, both atoms interact due to 
the cavity-stimulated exchange of a photon: |ei, a2, 0) 
|ai,a2,i) |<^i,e2,0) for a initially empty cavity |0). 
This interaction sequence contains in its middle part a 
virtual state and is independent of the photon number 
that was initially in the cavity. By following similar lines, 
therefore, the (initial) atomic state |ei, ^2, . . . , aAr) of 
atoms in the chain evolves according to the sequence 

|ei, a2, . . . , ttAT, 0) ^ 1^1, 62, ^3, . . . , «iv, 0) 
^ |ai,a2,e3, . . . ,aAr,0) 

|ai,a2,a3, . . . ,eAr,0), 

and where the virtual state |ai, . . . , aAr, 1) has been omit- 
ted for brevity. As seen from the sequence (|3]) , the atoms 
in the chain interact due to the cavity stimulated ex- 
change of a single photon between the originally excited 
atom and one of the N — 1 other atoms. This photon 
exchange, moreover, requires a rather large detuning be- 
tween the transition frequency of the atoms and the res- 
onant frequency of the cavity mode 

I {uje - ^a) - I > g{ri). i = 1, . . . , , (4) 

namely such that the cavity remains almost unpopulated 
in the course of interaction [16]. In the expression above, 
we have introduced the position-dependent atom-cavity 
couphng 

g{r) = go exp [-z^/w^] , (5) 

which is caused by the variation of the transversal cavity 
field along the atomic trajectories, and where go denotes 
the vacuum Rabi frequency and w the cavity mode waist, 
i.e. one half of the minimum width of the cavity field [see 
Fig.E^a)]. 

Recall that according to our scheme, however, the 
atoms are loaded to the cavity in the product state 
|li, a2, . . . , ttAr) and hence an intermediate excitation 
|1) \e) is first needed to bring the atoms to inter- 
action by means of the detuned cavity (see above). In 
order to realize this excitation, the atomic chain is ex- 
posed to a detuned laser beam with frequency ul that 
couples — via a position-independent strength ft to the 
|1) ^ |e) transition as shown in Fig. [2l With the cou- 
plings of the atoms to the both laser and cavity fields, 
the initial atomic state |li, a2, . . . , aAr) evolves according 
to the sequence of intermediate atom-cavity states 
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|li,a2, . . . ,aAr,0) |ei, a2, . . . , aAr, 0) |ai, a2, . . . , aAr, 1) |ai, 62, as . . . , aAr, 0) |ai, I2, as, . . . , aAr, 0) 

^ |ai,a2,e3,...,aAr,0) ^ |ai, a2, 13, • • • , ^at, 0) 

|ai,a2,a3,...,eAr,0) |ai, a2, as, . . . , Iat, 0) 



into the one of the final states |ai, I2, as, . . . , aAr), . • 
|ai, a2, as, . . . , Iat), which can have the same probability 
to occur. 

In Refs. [13I, it was shown that the condition 

\{lje-cji) -cjlI > ^ (7) 

for the detuning between the atomic transition and 
laser frequencies, ensures that the states |ei, a2, . . . , aAr), 
|ai, 62, . . . , aAr), |ai, a2, . . . , gat) remain almost un- 
populated. For this reason, the condition ([7]) plays the 
same role as the condition Q for the atom-cavity inter- 
action which makes the state |ai, . . . , aAr, 1) to be only 
virtually populated. In the following, we shall omit these 
unpopulated (intermediate) states and express the se- 
quence (j6j) in the short form 

|li,a2, . . . ,aAr,0) |ai, I2, as, . . . , aAr, 0) 
^ |ai,a2. Is, . . . ,aAr,0) 

|ai, a2, as, ... , Iat, 0). 

With the evolution (j8j) of the atomic chain due to the 
laser-cavity mediated interaction, the entangled W state 

N terms 

1 . , / s 

-7={e''^ |li, a2, . . . , aAr) + . . . + |ai, a2, . . . , Iat)) (9) 
V A' 

can be, in principle, generated by adjusting properly the 
atomic velocity v and the inter-atomic distance d for a 
given set of fixed cavity and laser parameters. In the 
last step of our scheme, the electronic population of each 
atom that leaves the cavity is (coherently) transferred 
back from the state |a) to the hyperfine state |0) in order 
to protect the atom against the spontaneous decay. This 
back transfer is achieved by applying the Raman pulses 
L2 behind the cavity and it produces the state ^ from 
([9]), once all the atoms have left the cavity. By utilizing 
the proposed scheme, therefore, a multipartite W state 
can be generated starting from the initial product state 
([2]) that is associated with the chain of N atoms which 
are conveyed through the cavity. 

In the next section, we shall analyze in details of how 
the laser-cavity mediated evolution (|8|) depends on the 
velocity v and the inter-atomic distance d when the 
atomic chain is conveyed through the cavity. In order to 
take into account these two parameters, we shall consider 
the position-dependent coupling (|5|), which gives rise to 



the time-dependent coupling between the i-th atom and 
the cavity 

g,{t)=goexp[-{z^^vt)yw^], i = l,...,N, (10) 

where and v denote the initial position of the i—th 
atom and its velocity along the 2;-axis, respectively. 

III. EFFECTIVE HAMILTONIAN AND 
MULTIPARTITE DYNAMICS 

While sequence (j6j) displays the basic concept of how 
the cavity-laser mediated interaction is achieved be- 
tween the atoms, we still have to analyze this cou- 
pling quantitatively as to understand how to control 
it in practice. For this purpose, we shall adiabati- 
cally eliminate the intermediate states |ai, . . . , aAr, 1) and 
|ei, . . . , aAT, 0), . . . , |ai, . . . , eAr, 0) from the sequence (|6]). 
This shall lead to an effective Hamiltonian that describes 
the time evolution of N atoms which interact with each 
other according to the simplified sequence (|8]). 

To outline this elimination process, let us first intro- 
duce the short-hand notation 

|Vi) ^ |V^+i) ^ |Vo) ^ |V^+2) ^ IV2) 

^ |V^+s) ^ |Vs) ^^^^ 

|V2iv) ^ |V^), 

for the composite states of TV identical atoms and the 
cavity, which corresponds one-to-one to the states from 
sequence (|6|). With this notation, the W state (|9]) refers 
to the states |Vi), . . . , IVat), while the cavity- mediated 
photon exchange is performed between the state IVat+i) 
and (one of) the states |VAr+2), • • • , |V2Ar), respectively. 

For N identical atoms, the evolution of the coupled 
atoms-cavity-laser system is described by the Hamilto- 
nian 

N 

H = uc c^^{ui\l)i{l\ ^ UE\e)i{e\ ^ UA\a)i{a\ 

+ [ne-^^-'\e),{l\^g,{t)c\e),{a\^h.c.]) , (12) 

where the first term describes the cavity energy, with 
c and c+ being the annihilation and creation operators 
for a cavity photon acting upon the Fock states |n), and 
(the summation of) the second term describes the chain 
of atoms and their interaction with the laser and cavity. 
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In this Hamiltonian, the interaction of the i— th atom 
and the cavity is based on the time-dependent couphng 
([To]). In the summation of the second term, moreover, 
each term contains the excitation energies cji, uoe^ and 
uja which correspond to the atomic states |e), and 
I a), respectively. 

In order to simphfy the analysis of the time evolution 
associated with the Hamiltonian (p!2]) . we switch to the 
interaction picture by using the unitary transformation 



Ui = exp ( -i^ (cji + cjL)t ^ |e)i(e| - Zcjit ^ 

\ z=l i=l 

/AT \ 

X exp ( —Luja \a)i{a\ — I [ujl — oja + ic^ c\ 



In this picture, the Hamiltonian (p!2]) takes the form 

AT 



Ri = -5c-^c^"^Hi with 

i=l 

Hi = A\e)i{e\ + [Q \e)i{l\ + gi{t) c \e)i{a\ + h.c] 



{13) 



and where A = {uje — (^i)—(^l and S = {ul—ojc) — {^a — 
uji) refer to the two off-resonance shifts (detunings) of the 
laser and cavity frequencies, respectively, as depicted in 
Fig. [2jc). The time evolution of the wavefunction then 
follows the Schrodinger equation 



(14) 



where we restrict the wavefunction |^) to the ansatz 

1^) = exp ( Z—t) ^a(t)|V,), C,(0) = fo. (15) 

For this ansatz, the Schrodinger equation (p!^ gives rise 
to a closed system of 2N + 1 equations (i = 1, . . . , A^) 



ZCo(t) 



^)^o(t) + ^^,(t)Civ+,(t), 



^Ci{t) = — a(t) + QCiv+i(t), 



Q2 



(16) 



LCN+i{t) = \^—+AJCN+i{t) + nCi{t) + gi{t)Co{t), 

and where the dot denotes the time derivative. 

As explained above, the + 1 states |Vo) and 
IVat+i), . . . , |V2Ar) remain (almost) unpopulated if the 
atom-cavity and atom-laser detuning satisfy the two con- 
ditions dlj and (O, respectively. Therefore, in order to 
separate these states from Eqs. (p!6|) . we utilize the adia- 
batic elimination procedure which assumes an adiabatic 
behavior of the functions Co{t) and CAr+i(t), . . . , C2N{t) 
or, in other words, that their time derivative vanishes to 
a good approximation. Together with condition ([7]) and 
conditions 

\5A\:^n\ |5A|»5(t)2, (17) 



it is justified to eliminate N-\-l equations from the system 
(p!6|) . Here, we shall omit further details of this derivation 
for which the reader is refereed to the seminal paper [l6| . 
The remaining N (effective) equations for the functions 
Ci(t), take the form 



N 



J = l 

(J#i) 



(18) 



Owning to the Eqs. (p!8|) for functions Ci{t) and the 
adiabatic behavior for functions Co{t) and CAr+i(t), . . . , 
C'2Ar(t), the evolution of the overall state of the atomic 
chain is given by the wavefunction 



i=i ^ -^-^ ^ 

which is associated with the effective Hamiltonian 

N 



(19) 



nit) 9 J it) 

26 



{StS-+S-S+), (20) 



where 5*^^ = |l)i(a| and 5*" = |a)i(l| denote the atomic 
two-photon excitation and de-excitation operators. Ob- 
viously, this Hamiltonian (|2Q]) is much simpler and de- 
scribes the effective atomic evolution (|8|) as mediated by 
the combined laser and cavity fields in Eqs. (fT8|) . In order 
to summarize the steps before, we have therefore found 
that the evolution of four-level atoms is reduced to the 
evolution of effectively two-level atoms which interact via 
a two-photon exchange in such a manner, that the excited 
state |e) remains (almost) unpopulated. 

When all the N atoms have left the cavity, the wave- 
function (p!9|) becomes in the limit t +oo 



/ + 00 
H,eds, (21) 
-OO 

where the matrix elements Mij = (Vi|M|Vj) are given 

by 

Mii=0 and Mi, =0{v,\i-j\d) for i j, (22) 
with 



0{v,d) 



TT gjw 
2 5A^v 



exp 



d^ 



(23) 



The latter expression ()23|) can be interpreted as the 
asymptotic coupling for a pair of atoms that move with 
the same velocity v and are separated by the distance d 
from each other. The atomic evolution of the state ([2T]) . 
therefore, is completely characterized by the asymptotic 
coupling (|23]) which depends on the two parameters (^;, d) 
once the frequency shifts and couphng constants ^, A, 
^o, and Vt are fixed by a particular experimental setup. 

Let us recall here that, after the atoms have left the 
cavity, their population in the (electronic) state |a) is 
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dist. units 



FIG. 3: (Color online) (a) Atomic velocities v and inter-atomic distances d for which the initial product state of two atoms 
1 11,02) becomes maximally entangled (solid lines). The dashed lines, in contrast, indicate the {v,d) pairs for which the atomic 
qubits remain disentangled, (b) Plot of the von Neumann entropy E{v,d) as a function of the atomic velocity and distance. 
In all these figures, the velocities v are displayed in units of and the distances d in units of w. 



transferred coherently back into the hyperfine state |0) 
by applying the Raman pulses L2. With this transfer, 
the wavefunction ([2T]) then gives the entangled W-class 
state of N hyperfine qubits 

N 

l*W) = E^^(^'^)l^i)' ^ith (24) 

i=l 

|Vi) = |li,...,0^), |V^) = |Oi,...,liv), 

and where the functions 

Ci{+^) = Ci{v,d) = (V:|exp(-tM) IV'i) (25) 

are obtained from the exponentiation of the hermitian 
operator (— ZM). These functions can be computed rou- 
tinely for any number of atoms for instance, by diag- 
onalization of the matrix (j22j). 

The wavefunction (|24|) , however, has not yet the de- 
sired form of a W state ([1]). In the next subsections, we 
shall therefore discuss the properties of the Wat states 
for different values of N and display those v and d pa- 
rameters, for which the function |^^) is equivalent (or 
close) to the desired W states. 

A. Two-partite entangled state 

For a chain of just two atoms (N = 2), the wavefunc- 
tion (|24|) takes the simple form [13| 

1^2) = cos 0{v,d) \\[) - L s'me{v,d) IV2) . (26) 

From this expression, we readily recognize that the two- 
partite W2 states 

|W^±) = i=(e±^f|li,02) + |0i,l2)) (27) 



are obtained (up to a global phase factor) if the condi- 
tion Oiv^d) = (2n + l)7r/4 is satisfied for some integer 
n. For a fixed set of experimental parameters (5, A, 
Qo^ and Vt^ therefore, a maximum entanglement is ob- 
tained only along the solid lines displayed in Fig. [H^a) for 
n = 0, 1, 2, 3. Obviously, the change between maximally 
entangled (solid lines) and completely disentangled states 
(dashed lines) happens more and more rapidly as the ve- 
locity of the chain is decreased from a certain maximum 
value (namely, for n = 0) onwards. 

Apart from understanding the dynamical parameters 
iv^d) for which a maximum entanglement is achieved, 
it is important also to know how sensitive these states 
are with regard to small uncertainties in the velocity 
and inter-atomic distance. To analyze this sensitivity. 
Fig. [3fb) displays the von Neumann entropy [18j 

E{v,d) = -Tr(plog2/9) 

= — cos^ ^(^;, d) log2 [cos^ ^(^;, d)] 

— sin e{v,d) log2 [s\n^e{v,d)\ ,(28) 

for velocities and distances satisfying 0{v^d) < 2 7r, and 
where p = Tt2 (1^2) (^2 1) denotes the reduced density op- 
erator with regard to the second hyperfine qubit. As ex- 
pected, the maximal values of the von Neumann entropy, 
i.e., E{v^d) = 1, are obtained along the lines which are 
displayed in Fig. [3l^a). Moreover, the least rapid varia- 
tion in the maxima occurs along the n = line and for 
rather small inter-atomic distances. For small velocities 
or some larger distance of the atoms, in contrast, a good 
control of the entanglements of the |<l>2) states becomes 
more and more difficult. 

Fig. [3fa) shows that an entanglement between the 
atoms occurs even for inter-atomic distances which are 
larger than 2w^ i.e. twice the cavity waist. In a high 
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finesse cavity, indeed, the Gaussian profile (j5j) approxi- 
mates quite well the intra-cavity field and, thus, it is pos- 
sible to generate an entangled state even for the atomic 
separation d > 2w. In practice, however, the cavity relax- 
ation and the spontaneous decay of the atoms introduce 
certain limitations on the distance between the atoms, 
beyond which it is not possible to generate the entan- 
gled state (|27|) . In order to estimate this limitation, we 
consider the condition 



(29) 



which ensures that N atoms couple strongly to the cavity 
field and, therefore, implies the validity of the effective 
evolution (|2Q]) . Here, and 7 denote the cavity loss rate 
and the atomic decay rate, respectively. For N = 2^ 
the above condition bounds the atomic coordinate to the 
interval Z- < z < Z-^ with 



± w 



ln[2 g^/iKj)] 



(30) 



Owning to these boundaries, therefore, the distance d 
between two atoms must satisfy 



d 



-< V2 ln[2 52/(K7)] = ^- — 

W WW 



(31) 



For the typical atom-cavity parameters [l2|: {^o, 7} = 
27r X {10,0.4,2.6} MHz, this codition implies the lim- 
itation d < 3.243 w. We note that this estimation 
agrees well with the solid lines from Fig. [3fa) since, for 
d > 3.2 w^ the atomic velocity becomes so small that it 
would prevent any experimental implementation of our 
scheme. 



B. Tree-partite W state 

For a chain of three atoms {N = 3), the wavefunction 
241) takes the form 



m = j2^'.iv,d)\vi), 



(32) 



i=l 



with 



Ci{v,d) = 

C2{v,d) = 
C3{v,d) = 



6 : 



(33) 




1.0 1.5 
veloc. units 



2.0 





FIG. 4: (Color online) (a) Lines along which the moduli 
\Ci(v,d)\ (solid), \C2(v,d)\ (dashed), and \C3{v,d)\ (dotted) 
are equal to 1/V3. These lines correspond to velocities ()36p 
with n = 0, 1 and m = 1, 2 (in the limit d ^ 0). (b) Two max- 
ima in the fidelity F-{v,d) obtained for velocities (|36)) with 
n = 0, 1 and m = 1,2 (in the limit d ^ 0). For guidance of 
the eyes, the semi-transparent layer displays a constant value 
F-{v,d) = 0.5. (c) The same as in figure (b) but for the fi- 
delity F-^(v,d). Again, all velocities v are displayed in units 
of glw/5 and the distances d in units of w. 
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and where we used the notation ^ = exp [— (i^/(2 k;^)] 
and 



A4 



exp 



±1, 



8 SA^v ' 

c = ix{e + v»+e) ■ 

In order to obtain the state W3 from wavefunction 1^3), 
we have to determine those pairs (v^d) for which the 
equations 



\Ci{v,d)\ = \C2{v,d)\ = \Cs{v,d)\ 



1 

71 



(34) 



are fulfihed. In Fig.SJ^a), we displayed the corresponding 
Hues for which the moduh \Ci{v^d)\ (sohd), \C2{v^d)\ 
(dashed), and \Cs{v^d)\ (dotted) are equal to l/>/3. The 
requested W3 states are obtained for those (v^d) pairs, 
for which all three types of lines intersect with each other. 

As seen from Fig.m^a), however, the lines for the (mod- 
uli of the) functions Ciiy^d) intersect only if the inter- 
atomic distance vanishes. In order to determine the cor- 
responding velocities, we first observe that for — > 
(^ 1), the wavefunction (|32]) becomes 



^ 1 ^ 2 p*^^ 



z4x 



1 



I V2) + |v^)) 



3 ' 3 
This expression can be readily cast into the form 



W-3±) = i=(e±^^|V;) + |Vi) + |V^)) 



(35) 



if X = (3 n + m) 7r/9 or, equivalently, if the velocity takes 
the values 



8 5/S? 7r(3 n + m)' 



(36) 



with m = 1, 2 and n being an integer. To summarize, the 
vanishing inter-atomic distance along with velocities ([36|) 
are both necessary to obtain the W3 states due to wave- 
function (|32]) . According to our scheme, however, the 
atoms are separated by a macroscopic distance which is 
non- negligible with regard to the cavity waist. Therefore, 
we may determine the parameter region {v^d) with the 
non-zero inter-atomic distance, for which the two fideli- 
ties [13 



F±{v,d) 



mtw,)? (37) 

^ ^ Cl(^;, d) + C2{v, d) + C^{v, df 



between the states \W^) and l^g) are larger than the 
threshold value of 1/2. These fidelities F-(v^d) and 
Fj^(v^d) are displayed in Figs. (|l])(b) and (c), together 




FIG. 5: (Color online) Fidelity F4{v,d) for the generation of 
the W4 state as a function of the velocity and the inter-atomic 
distance. Again, the maximum value F4{v,d) = 1 is obtained 
only for vanishing distance (d = 0) and velocities (|4Qp with 
n = 0, 1. The semi-transparent plane with F^iv^d) = 0.5 is 
plotted to guide the eyes; the units are the same as in Figs. [3] 
andSl 



with a semitransparent plane in order to delimit the re- 
gions for which F±{v^d) > 0.5. While the maximum 
values F±{v^d) = 1 are obtained only for a few veloc- 
ities and vanishing inter- atomic distance, there are still 
{v^d) regions (with non-zero distance) for which the fi- 
delities become reasonably close to the maximal value. 
Note, moreover, that the region with Fj^{v^d) > 0.5 is 
notably larger than those with F-{v^d) > 0.5. We con- 
clude, therefore, that from the experimental perspective 
it might be preferable to generate the \W^) state be- 
tween three hyperfine qubits by means of the suggested 
scheme. 



C. Four-partite W state 



For a chain of four atoms {N 
dsn can be written as 



4), the wavefunction 



W^) = Y.C,{v,d)\V[). 



(38) 



In contrast to A/" = 2 or = 3, however, the expres- 
sions for Ci{v,d) become rather bulky now and are not 
displayed here. Recall from the previous subsection that 
the wavefunction |<I>3) produced the \W^) states only 
for vanishing distances and velocities ([36|) . In this sub- 
section, therefore, we proceed in a similar fashion and 
consider the wavefunction ([38|) with the vanishing inter- 
atomic distance {d 0) 



-i4x 



l + 3e^^^ 



VI) 



-i4x: 



^8x 



From this expression, the state 
1 



|w^4) = T.(e*nv;) + |v^) + |v;,) + |v;)) 



(39) 
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is readily produced if x = 7r(2 n 
if the velocity takes the values 



l)/8 or, equivalent ly, 



IT Vt^gf^w 8 

8 7r(2n + l)' 



(40) 



More generally, Fig. (|5]) displays the fidelity ^4(1^, d) = 
|(W4|$4)P for the generation of the W4 states as a func- 
tion of the velocity and the inter-atomic distance due to 
the wavefunction (|38|) . Analogue to the last subsection, 
the maximum fidelity F/^iv^d) = 1 is obtained only for 
zero distance {d = 0) and velocities that fulfill the condi- 
tion (|4Q]) . For non-zero distances, nevertheless, there is 
one broad parameters region for which the W4 state (|39|) 
can be generated with a reasonable hight fidelity. 



D. N >5 partite W state 

For any other number > 5 of atoms in the chain, the 
functions Ci{v, d), . . . ^ Cn{v^ d) can still be computed 
from the formula (|25]) . and the wavefunction can be 
further analyzed with regard to (v^d) region, for which 
the corresponding Wn state is produced most reliably. In 
order to generate such state, according to the definition 
([1]), the conditions 



\Ci{v,d)\ = ... = \CN{v,d)\ = 



1 



N 



(41) 



need to be fulfilled. By performing numerical analysis, 
however, it turns out that the Eqs. (|4T]) cannot be fulfilled 
for any choice of the velocity and inter-atomic distance 
and, therefore, the fidelity |(WAr|^^)P is always smaller 
than unity. 

Nevertheless, we can display this fidelity as a function 
of V and d and determine the region where it takes the 
highest value. In order to proceed so, however, we still 
need to specify the reference state IVKat), which we are 
looking in the form 



\Wn) 



1 



N 



N \ 

Vl> + ElV:> (42) 

i=2 J 



with an unknown phase (j). The form of this state has 
been chosen in line with the previously obtained W states 
([27|) . (j35j), and (|39|) . In order to calculate the unknown 
phase (j) in (|^2]) . let us first consider the W-class state 



N 



(43) 



N 



^ l + (4iA'-l)exp(Z 2Nx) ,,,, . 



N 



where the velocity v is such that it makes the expressions 



|C.KO)|--^ 



\Cn{v,0)\-^ 



(44) 



(a) 



F(N) 



10 15 20 25 30 

N 




FIG. 6: (Color online) (a) The fidelity F{N) = \{Wn\^n)\'^ 
has its maximum value F(N) = 1 for N = 4 and decreases 
monotonically as the number of atoms increases. For A > 15, 
it falls below the threshold F{N) = 1/2 (dotted line), (b) 
Fidelity F^{v^d) for the production of the \W^) state due to 
1^5) as a function of the velocity and the interatomic distance. 
Similar as in Fig. O it reaches its maxima for d — and 
velocities (|45|) with n — 0, 1. The units are the same as in 
Figs. [3] and |4] above. 



minimal. It can be straightforwardly shown that all the 
expressions (|44|) are minimized for the values x — ^ + 
l)/(2 A/") or, equivalently, for velocities 



8 (2n + 1) 



(45) 



Substituting this value for x, the state ([43|) then becomes 



N 



N 



N 



(46) 



For a vanishing inter-atomic distance, therefore, the state 
with velocities (|45)) gives the best approximation to 
the Wat state. 

As we explained in the beginning of this subsection, 
there are no such iv^d) pairs for which the Eqs. (|^T]) 
can be fulfilled. However, we found the state (|46|) which 
gives the best approximation to the Wat state (|^ . By 
comparing the states \Wn) with for A' = 4, we 
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find that the phase (j) is equal to tt and, therefore, the 
reference state becomes 

\WM) = ^(^e^''\V[) + f2\V'^y N>4. (47) 

In order to understand how well the state (|^6|) approx- 
imates the above state (|47|) , in Fig. [6ja) we displayed 
the fidelity F{N) = \{Wn\^n)\'^^ As seen from this fig- 
ure, the fidelity has its maximum value F{N) = 1 for 
N = 4 and decreases monotonically as the number of 
atoms is increases in the chain. The fidelity drops below 
the threshold F{N) = 1/2 for AT > 15. We therefore con- 
clude that the state (|^ approximates reasonably well 
the reference state for at least 4 < < 15. 

Having specified the reference state (|T7|) . we can eval- 
uate the fidelities 

FN{v,d) = \{WN\<i>'N)\^; 5<N<15 (48) 

as functions of the velocity and the inter-atomic distance. 
In Fig. ini^b), for instance, we display the fidelity (|48|) 
for N = 5. According to this figure, moreover, the fi- 
delity reaches its maxima F^iv^d) = F(b) ~ 0.97 for 
(i = and velocities that satisfy the condition (|^ with 
n = 0, 1. Let us note here that the typical spacing be- 
tween two neighbored potentials wells (sites) of an optical 
lattice is in the sub- micrometer range [1O-I2]. As seen 
from Fig. [6l^b), this typical spacing is comparable to the 
inter-atomic distance for which the fidelity F^{v^ d) ^ 0.9 
is reasonable high and where the typical cavity waist 
{w = 20 /im) has been considered as the distance units. 
The recent developments in cavity QED, therefore, make 
it possible to generate the W5 state by means of the 
proposed scheme. If we compare, however, the (v^d) 
regions for which the fidelities F±{v^d) [Fig. |4l^b)-(c)], 
F4{v,d) [Fig. [5], and F^{v^d) [Fig. [6l^b)] are higher than 
the threshold value of 1/2, we also conclude that these 
regions become smaller as the number of atoms (in the 
chain) increases. 

IV. REMARKS ON THE IMPLEMENTATION 
OF THE PROPOSED SCHEME 

In our discussions so far, we have always assumed 
that the velocity and the distance of the atoms in the 
chain, i.e. their position within the optical lattice, can be 
controlled exactly. With this assumption in mind, the 
atom-cavity coupHng was described by expression (p!Q|) . 
This assumption, however, neglects the transversal com- 
ponents of the cavity field as well as the oscillations of the 
atoms within the potential well due to their finite tem- 
peratures, which include both the axial (along the z-axis) 
and radial (along the x, ?/-axes) oscillations [see Fig.[2ja)]. 
This additional motion gives rise to a dispersion of the 
atomic positions and velocities and, thus, leads to un- 
certainties in selecting the dynamical parameters in our 
model. 



Obviously, any significant uncertainty in the parame- 
ters {v^ 5^ A, ^o, Vt} will influence the generation of the 
desired W entangled states. According to our scheme, 
however, these entangled states are produced when all 
the atoms have left the cavity. Instead of understanding 
these parameters as exacts therefore, they should refer 
to the mean values and we need to analyze how small 
(but realistic) variations in these parameters affect the 
final state of the atomic chain. For instance, the radial 
oscillations of the atoms lead to the mean value of the 
vacuum Rabi frequency 'g^ and axial oscillations to the 
mean values of the inter-atomic distance d and velocity 
U, respectively. Axial oscillations affects also the initial 
position z"- and velocity vi of each atom inside the lat- 
tice and result in uncertainties l\di = d — — and 
AVi = V — Vi, where i = 1, . . . , A^. Therefore, the plots 
E{v,d), F±{y,d), F^iv^d), and F^{v,d) from Figs. (|3j)- 
should be re-calculated as function of the mean values 
V and d and their corresponding uncertainties, respec- 
tively. 

In order to determine realistic uncertainties for the dis- 
tance and velocity of the atoms in the chain, we first 
mention that recent cavity QED experiments allow to 
position the atoms relative to the cavity antinode with 
a precision of ~ 0.1 /im by utilizing an additional dipole 
trap acting along the cavity axis [l^. When compared 
to the typical cavity wavelength (~ 0.8 /im), such posi- 
tioning precision leads to the spatial dispersion which, 
in turn, yields the mean value 'g^ ^ 0.75 go that is still 
good enough for our purposes [see Eq. (p!7|) ]. Moreover, 
the same spatial dispersion implies upper bounds for the 
uncertainties \Ad/d\ < 0.2 and \Av/v\ < 0.1, if compared 
with the typical spacing (~ 0.5 /im) between two neigh- 
bored potential wells of an optical lattice and the typical 
atomic velocities 0.5 m/s) along the lattice axis. 

For a further analysis of how reliably a given (experi- 
mental) setup will generate a particular W state, in Fig. [71 
we display the (mean) functions E{v^d)^ F±{v^d), and 
F4^{v^ d) by calculating their average for a certain spread 
of parameters. For each subfigure [71(a)- (d), we have ran- 
domly chosen 20 uncertainties Ad and Av from the inter- 
vals [0, 0.2 (i] and [0, 0.1 U], respectively. By comparing 
the Figs. [3jc) and[7fa) it can be seen that the von Neu- 
mann entropy, for instance, is reduced considerably for 
its sharp maxima (n = 3) and that it remains almost 
the same around the broad maxima (n = 0). Similarly, 
the mean fidelities which are displayed in Figs. [71(b)- (d), 
are considerably reduced for their sharp maxima. These 
(^;, d) regions for the velocity and inter- atomic distance in 
the atomic chain are, therefore, less useful for any prac- 
tical implementation and only the {v^d) regions which 
correspond to the broad maxima of the von Neumann 
entropy and fidelities, are relevant for the generation of 
entangled W states by means of the proposed scheme. 

Finally, we assumed in our treatment that the (center- 
of-mass) position of each atom is not a quantum variable 
but described classically by the vector fi{t) = {0, 0, + 
vt}. Obviously, such an assumption excludes several im- 
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FIG. 7: (Color online) Entanglement and fidelity measures averaged over 20 randomly chosen uncertainties Ad and Av of the 
inter-atomic distance and the velocity, respectively, (a) Von Neuman entropy E{v,d), (b) fidelities F-(v,d) and (c) F-^(v,d) 
for a chain of three atoms, and (d) F/^{v^d) for a chain of four atoms. The uncertainties for the distance and velocities are 
chosen from the intervals [0, 0.2 d] and [0, 0.1 iJ], respectively (see text for further details). 



portant effects on the atomic motion that arise due to 
quantization of the cavity field. For example, the corre- 
lations between the internal dynamics of the atoms and 
their transverse (center-of-mass) position may lead to an 
additional source of decoherence and disentanglement in 
the effective evolution (|2Q]) . if the atom-cavity system is 
embedded in a realistic reservoir [19]. In our scheme how- 
ever, the external potential that is created by the optical 
lattice dominates the kinetic energy associated with the 
atomic momentum. Therefore, any correlations which 
are induced by the mechanical effects of cavity on atoms 
are strongly suppressed and have been neglected in our 
treatment. 



V. SUMMARY AND OUTLOOK 

A scheme is proposed to generate the entangled W 
states for a chain of N four-level atoms that are equally 
separated and conveyed through an optical cavity by 
means of an optical lattice. This scheme is based on 



the cavity-laser mediated interaction between the atoms 
which are separated by a macroscopic distance, and 
works in a completely deterministic way for qubits en- 
coded by two hyperfine levels of the atoms. Only two 
parameters, namely the velocity v of the chain and the 
inter-atomic distance d, determine the effective inter- 
action among the atoms and, thus, the degree of en- 
tanglement that is obtained for the overall chain. The 
asymptotic coupling (|23]) that completely characterizes 
the atomic evolution, tells explicitly how the degree of 
entanglement depends on these two parameters. The 
purpose of this work is to understand the state evolu- 
tion of the atomic chain and how it can be utilized to 
generate the entangled W states. For chains consisting 
of A/" = 2, 3,4 and 5 atoms. Figs. ([3])-(|6]) display the von 
Neumann entropy and the fidelities as functions of the 
velocity and inter-atomic distance. For 5 < TV < 15, 
moreover, we suggested the reference state (|47|) which 
is approximated by the wavefunction |^^) with a high 
fidelity. In view of the recent developments in cavity 
QED, moreover, we have also analyzed and discussed the 
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proposed scheme with regard to sensitivity in the forma- 
tion of desired entanglement due to uncertainties in the 
atomic motion. 

For two or more atoms, the generation of entanglement 
by means of a (detuned) optical cavity has been investi- 
gated in several papers [ij, [H, [2Q|. All these studies, 
however, relied on the small sample approximation in 
which the separation of the atoms is considered to be 
negligible when compared with the cavity waist. Only 
recently (2ll-[23j. the atom-cavity coupling (|5]) has been 
exploited in more detail in order to suggest various en- 
tanglement schemes within cavity QED. In the work by 
Amniat-Talab et al., for instance, a scheme was proposed 
in which two atoms were coupled sequentially to a reso- 
nant cavity and where a position-dependent coupling is 
used to drive a STIRAP-type process in order to reduce 
the losses due to atomic and cavity decays. Moreover, the 
scheme by Marr et al. is also based on a STIRAP-type 
process and describes an adiabatic evolution of a prod- 
uct state of two atoms which are coupled simultaneously 
to a detuned cavity. The success of this scheme, how- 
ever, relies strongly on the ability to detect the photons 
which leak through the cavity mirrors with an efficiency 
close to one. In both schemes, therefore, the atomic ve- 
locities and inter-atomic separation are used to control 
the accuracy of a STIRAP-type process, in contrast to 
our approach, in which these parameters are utilized to 
control the degree of entanglement. 

Our proposed scheme might be suitable also for ion- 



cavity experiments in which N trapped ions interact si- 
multaneously with a (detuned) optical cavity [25| . In 
these experiments, the same coupling to the laser and 
cavity fields applies for ions with a three-level A-type 
configuration [as displayed in Fig. [2fb)]. For such a level 
configuration, the qubit is associated with the states |1) 
and I a), and the W state can be generated by moving the 
equally distanced trapped ions along the trap. Similar as 
for the atomic chains above, the cavity-laser mediated 
interaction between the ions is described by the effec- 
tive Hamiltonian (|2Q]) and, therefore, requires the same 
analysis as performed in Sects. III. A - III.D in order to 
produce the Wat states. 

Finally, we like to mention our recent paper [26] in 
which another deterministic scheme for generation of 
the multipartite W states has been proposed. In con- 
trast to this work, the resonant atom-cavity interaction 
regime has been exploited for N flying two-level (Ryd- 
berg) atoms which couple sequentially (one after another) 
to the mode(s) of a high- finesse bimodal cavity. 
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